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Abstract 

Nonfactorizable contributions to exclusive two-body nonleptonic weak decays of heavy 
mesons arising from color octet currents, characterized by the parameter r2, are extracted 
from the data. It is found that r2 is equal to ~ —0.67, —(0.9 ~ 1.1), — (1.2 ~ 1.3), 
~ 0.36 respectively for D —>■ Ktt, K*tt, K*p, B —>■ Dtt decays. As expected, soft-gluon 
effects become stronger when the decay particles are less energetic, allowing more time for 
significant final-state strong interactions. As a consequence, the parameter 02 is not universal 
and is channel or class dependent, contrary to the anticipation of the factorization approach. 
The leading 1/Nc expansion works most successfully for D Ktt decays as the sub leading 
1/Nc factorizable contribution is almost compensated by the soft-gluon effect. We argue that, 
in contrast to what happens in ^ D^*^Tx{p) decays, the nonfactorizable term and the 
subleading l/Nf. factorizable term in 5 — ipK'^*^ decays are opposite in signs, in accordance 
with a recent QCD sum rule calculation. Implications are discussed. 
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1. Introduction A common wisdom for having an approximate description of the under- 
lying physics of nonperturbative QCD is to search for the small parameters in the theory, 
then set the small parameters to zero and treat the finite effects of the parameters as small 
perturbations [1,2]. Hence, even for a nonperturbative physics such as QCD, one can still 
perform a sensible perturbative calculation. The well-known examples are chiral perturba- 
tion theory and the heavy quark effective theory which respresent a good approximation 
of QCD in the light and heavy quark mass limits, respectively. The corresponding small 
expansion parameters in these theories are ^g/AgcD and Aqcb/i^q with (mg) being the 
light (heavy) quark mass. 

By treating the number of color degrees of freedom as a free parameter Nc, 1/Nc can be 
another useful small parameter in QCD. The fact that l/N^ is indepedent of energy has its 
advantage and disadvantage. On the one hand, it can be used as an expansion parameter 
for both short- and long-distance dynamics owing to its energy independence. On the other 
hand, whether or not the 1/Nc approach works is at best case by case dependent, since, 
unlike chiral perturbation theory or the heavy quark effective theory, there is no certain 
kinematic region where the validity of the l/N^ expansion is guaranteed. Empirically, while 
the large-A'^c approach fails in some cases, it explains well the OZI rule and even provides 
quantitative predictions on the relative strength of the baryon-meson coupling constants 
and on the baryon mass relations up to 0{1/N^), as was realized last year [3]. It is also 
known for sometime that the leading 1/Nc expansion operates reasonably well for exclusive 
two-body nonleptonic decays of charmed mesons [4] ; the discrepancy between theory and 
experiment for color-suppressed channels e.g. —>■ K^tc^ is greatly improved provided 
that contributions from Fierz-transformed currents, which are suppressed by order l/N^ are 
dropped [6,7]. Due to the success of the l/Nc approach to charmed meson decays, it has 
been widely believed by many practitioners in this field that it applies equally well to the 
weak decays of bottom mesons. 

The recent CLEO data on the B decays B Dn, Dp, D*7t, D*p exhibit a quite striking 
result [8]: The interference between the two different amplitudes contributing to exclusive 
two-body B^ decays are evidently constructive, contrary to what naively expected from 
the leading 1/Nc expansion. In other words, the observed destructive interference pattern 
^For a review of the expansion for nonleptonic weak decays of mesons, see Ref.[5]. 
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in K^TT^, K^P^, K*^7T^ cannot be extrapolated to B decays. Though, as stressed 

before, the rule of retaining only the leading term in the 1/Nc expansion is case by case 
dependent, an understanding of why it is operative for charm decays but not for the B 
meson case is called for. 

Since A^^^ = 3 in reality, when the 1/Nc approach works empirically, it will imply either 
that the effective expansion parameter is something like l/(47rA'"c), 1/N^..., or that there 
is a dynamic reason for the suppression of non- leading l/N^ terms. In the large- A^^^ limit, 
the meson decay amplitude is factorizable [4] (this is no longer true for baryon decay as the 
baryon contains Nc quarks.) Irrespective of the 1/Nc expansion, the factorization hypothesis 
means that the meson two-body decay amplitude may be expressed as the product of two 
independent hadronic currents. Consider the operator 

Oi = {sc){ud), (1) 

where (^1^2) = (^ll^i{l — 15)12 ■ The decay amplitude of D'^ — > K^tt~^ induced by Oi in general 
can be expressed in terms of factorizable and nonfactorizable contributions: 

{K\+\0,\D^) = {n+\{ud)\0){K'\{sc)\D+) + {K\+\0,\D+)^f, (2) 

where the subscript nf denotes a nonfactorizable contribution. In the 1/Nc expansion, the 
factorizable amplitude is of order N^^'^, while the nonfactorizable one is of order N'^/"^. 
Using the identity 

Oi= ^02 + \{-s\''d){u\'^c), (3) 
with O2 = {sd){uc) and (giA"g2) = lilnO- ~ 'l5)^"'Q2, one obtains 

{K\+\0,\D+)^f = j^{K\+\02\D+)f + ^{K%+\{sX^d){uX^c)\D+) 

+ ^^{K\^\0,\D^U, (4) 

where the factorizable amplitude {K^tt^\02\D^) f is equal to {K^\{sd)\0) {7t^\{uc)\D^) . In 
the traditional vacuum insertion method only the factorizable terms are retained so that 

{K%+\Or\D^U = ^^{K%^\0^\D+)f. (5) 
^For example, the expansion parameter in QED is a = e^/47r even though e ^ ^ is not very smalL 
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However, this naive method encounters two problems. First, it is not logically consistent [9]: 
A nonfactorizable term becomes a subleading 1/Nc factorizable contribution. Second, naive 
factorization cannot explain the bulk of the experimental data of charm decay (for a review, 
see Ref.[5]). Both problems hint that it is indispensible to take into account nonfactorizable 
effects. To the order 0{N~^l'^), we thus keep the second nonfactorizable term on the r.h.s. 
of Eq.(4) and drop the third term as it is further suppressed by a factor of To the 

end, we have 

{K\^\0,\D^) = {K\^\0,\D^), + {K\^\0,\D^), + f ) , (6) 

where 

_ (^07r+|(sA"d)(iZA"c)|i:>+) 

" {KH+\{sd){uc)\D+) f ' ^ ' 

first introduced in Ref.[10], denotes the nonperturbative effects arising from the soft gluon 
exchange between the color octet currents (sA^c?) and (-uA^c) relative to that from the cor- 
responding color singlet currents and is of order l/N^ in the large- A'c limit. 

It is clear from Eq.(6) that if the leading 1/Nc expansion works, it will come from the 
dynamical reason that the subleading 1/Nc factorizable contribution is compensated by the 
nonperturbative correction, namely r ~ —2/Nc. Though in practice it is very difficult to 
estimate the nonperturbative soft gluon effects (for some recent attempts; see Refs. [11-13]), 
they can be extracted from the available data. Consider the decays D — > PP, VP, VV 
and B PP (P: pseudoscalar meson, V: vector meson). It is expected that the soft-gluon 
effect is such that 

\r2{D ^ VV)\ > \r2{D ^ VP)\ > \r2{D ^ PP)| > |r2(P ^ PP)|, (8) 

as the final-state particles have more time to allow significant final-state strong interactions 
when they become less energetic. The purpose of this paper is to determine the parameter 
r2 from data and confirm the above expectation. Implications on the factorization method 
are discussed. 

2. Charmed Meson Decays In the large- A^c approach both soft-gluon nonperturbative 
effects and final-state interactions are subleading 1/Nc nonfactorizable corrections [4]. In 
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order to determine the soft-gluon effects, we will focus on the exotic channels e.g. 
K^7i~^, n^TT^ and the decay modes with one single isospin component, e.g. D'^ 7r^</), —>■ 
7r+0, where final-state interactions are presumably negligible. 

The QCD-corrected weak Hamiltonian for Cabibbo-allowed charm decays is given by 

Kfr= ^V:,VUciO, + C20,), (9) 

with Oi = {sc){ud), O2 = {sd){uc), Ci and C2 are Wilson coefficient functions determined at 
the renormalization scale fi ~ rric 

ci(mc)~1.26, C2(mc) 0.51. (10) 

The amplitude for the decay — > K^7t~^ up to the subleading order in 1/Nc is 

a= a,{ml-ml)Uf^^{ml), (11) 

where = 132 MeV, fx = 160 MeV, 

Ol = Ci + ^iC2, 02 = C2 + 6ci, (12) 
with ^1 = 1/Nc + ri/2, ^2 = 1/N, + r2/2, and 



ri{D Ktt) 



{7r+K%ud) {sc)\D+)f {71+ R- \ (ud) (sc) 



/ 



- {n+K^\{sX''d){uX''c)\D+) _ {7r^K°\{sX^d){uX''c)\D^) 

r2[D ^ Kit) - (^^+KO\^sd){uc)\D+) f ~ {7i0K0\{sd){uc)\D0) f ' ^ ' 

where we have applied Eq.(6) to derive (11) and isospin symmetry to relate the K^n^ 
matrix elements of color octet currents to that of K^n^, K^tt^. In the literature, 

factorization often means that the parameters ai and 02 (and hence ^1 and ^2) are universal, 
namely they are channel independent in D or i? decays. However, we see from Eq.(13) 
that a priori there is no reason to expect that ri and r2 are decay mode independent. 

For the form factor f^^{q^) in Eq.(ll) we will use the average value 

f^^iQ) = 0.76 ± 0.02 (14) 



■^Note that = = corresponds to naive factorization, while leading \/Nc expansion leads to 

6=6 = 0. 
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extracted from the recent measurements of D ^ Kiu by CLEO II, E687 and E691 [14]. 
As for the form factor f^'^, there are only two available experimental information. An 
earher measurement of the Cabibbo-suppressed decay L>° — > 7r~£~^u by Mark III yields 
f+''{0)/f+^{0)\ = l-Otoi ± 0.1 [14,15], while a very recent CLEO-II measurement of 
D+ n°e+iy gives |/:^^(0)//:^^(0)| = 1.29 ± 0.21 ± 0.11 [16]. Though the latter perfers 
a larger f^'^{0) over f^^{0), its error is still very large. Fortunately we can use the recent 
measurement of D+ 7r+7r° by CLEO II [17] to fix /:^'^(0). The amplitude of D+ 7r+7r° 
reads 



AiD^ tt^tt') = ^ W.(ai + a2){ml - mDU^^mH (15) 



where ai and 02 are defined in the same manner as in (12) except that in the present case 
ri{D tttt) = r2{D ^ tttt). Assuming a monopole behavior for the form factor 

- 1 ^2/ 2V (16) 

1 - [q^/mi) 

where m* is the mass of the low- lying 1~ resonance that couples to the weak current, and 
using the experimental value [17] 

B{D+ 7r+7r°) = (0.22 ± 0.05 ± 0.05)%, (17) 

as well as t{D+) = 10.66 x lO'^^^ [18], we find [19] 

[l + CiD^7r7r)]f^^iO)^0.83. (18) 

To proceed further we assume that ^i{D Kn) = C,2{D ^ Kn) (which is at least valid in 
the SU(3) limit)= ^{D tttt). Substituting (14), (16), (18) into (15) we find ^{D Kn) ^ 
7 X 10"^, and hence 

r{D Ktt) f=i -\ . (19) 

Evidently, the subleading 1/A^c factorizable contribution is almost compensated by the non- 
factorizable soft gluon effect, so that ^ ~ 0. This explains why the leading 1/A^c expansion 
operates for D Kn decays. 

We digress for a moment to make a remark on the IV-exchange amplitude in the decays 
— > K~7r~^, K^n^, which is given by 

A(iy - exchange) = ^V:,V,s{c2{K-7i+\{sd)\0){0\{uc)\D') 

+ Ci(X-7r+|(sA"ci)(iZA"c)|L'°)}. (20) 



We see that although the first term vanishes in the SU(3) hmit owing to the conservation 
of vector current, the effect of soft-gluon exchange could be important. Indeed, a general 
phenomenological analysis of D — > Kn data indicates that 1^-exchange is small compared 
to external and internal 1^-emission amplitudes but it is not negligible [20]. 

We next determine the parameters Vi and r2 for D — > VP decays, whose amplitudes are 
of the form 

A{D^VP)= M{e-pD): (21) 

where £^ is the polarization vector of the vector meson V, and po is the momentum of the 
charmed meson. We first consider the decay — > ^tt"*" which proceeds solely through 
internal W emission: 

Ap+-07r+) = ^V:,V,sa2mrssm{n^\{uc)\D^) 

= ^V:,V^M{e-pn), (22) 

with 

b' ^ 2a2Um^f^^{ml), (23) 

where we have applied the relation (0|(ss)|O) = f^m^e^. Prom the measured branching ratio 
B{D+ 07r+) = (6.0 ± 0.8) x lO'^ [18], and the decay rate formula 



r(D+ ^ 07r+) 

we obtain 



' _ „A , (24) 



\h\D+ 07r+)| ^ 0.336 GeV^ (25) 

Comparing (25) with (23) leads to (/</, = 0.221) 

^2(^+ ^ 07r+) ^ -0.12 , r2{D+ 07r+) ^ -0.92 . (26) 

It is evident that the soft-gluon effect is larger than that in D — >■ PP decay owing to the 
small relative momentum between the VP final states. 
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The decay K*^7i^ receives both external and internal VT-emission contributions: 

A{D^ K*\+) = ^V:,V,s {a,{n^\iud)\0){K*\7sc)\D^) + a,{K*'\i-sd)\0){n^\iuc)\D^) 
= ^V:,VUa' + b'){e-pn), (27) 

with 

a' = aif^[{mD + mK*)A^^\ml) - {niD - mK')A^^* {ml)], 

h' = 2a2fK'mK*f+Hml,), (28) 

where use has been made of 

{V{pv)\{sc)\D{pd)) = i{e,{mD + mv)A,{q')-e-q^^^^-^A2{q') (29) 

mD + my 

- 2'-^q,mv[A,{q') - Ao{q% + ^ e^,„^e^p>^\/(g2), 



with q^ = {pd — pv)fj. and ^3(0) = Aq{0). Assuming ,^1 ^ ,^2 and using the experimental 
branching ratio B{D^ K*^7t^) = (1.9 ± 0.7)% [18], and the measured form factors [14] 

^' (0) = 0.59 ± 0.04 , A^^' (0) = 0.44 ± 0.09 , (30) 

and fx* = 0.220 determined from the decay r K*"^^-, we get Q 

^(D ^ ir*7r) ^ -0.22 , r{D ^ K'-tt) ^ -1.10 . (31) 

Thus far we have assumed .^1 = .^2 or ri = r2. Whether this is true or not can be tested 
from the decay —>■ 07r^ which occurs solely through the external W^-emission diagram 

A{Dt - </)7r+) = 9Zv:,V,,a'{e-pJ, (32) 

with 

a' = aMmo^ + m^)Af^*{ml) - {mo^ - m^)A^^'^{ml)]. (33) 
From B{D+ 07r+) = 3.6% [21] and r{D+) = 4.50 x lO^^^^ [18], we find 



\a'\ = 0.25 GeV^ (34) 
^For some unknown reasons, the result for ^ determined from the decay Z)+ p'^ is very unsatisfac- 
tory: ^{D Kp) 1.4. 
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If the form factors taken to be the same as Af2* : we then obtain 

^i{D, (f)Tr) ^ -0.27 , ri{Ds (pn) ^ -1.22 . (35) 

Since ^1,2 = (01,2 — Ci,2)/c2,i and \c2\ << \ci\, it is clear that the determination of ^2 is far 
more less uncertain than ^1. Within uncertainties we see that ri and r2 are very close. Note 
that internal W emission is usually suppressed relative to external W emission owing to color 
mismatch. However, we see from (25) and (34) that internal W emission in D — > 07r, D — > 
Ktt dominates over external W emission. Of course, nonfactorizable soft-gluon contribution 
plays an essential role here. 

For the decay D — > VV, its general amplitude reads 

A{D ^ V1V2) = i^V:^Vg^,^e^iV^)e''iV2)iA,g,, + A^p^^p^ + tVe,,apP%P^i). (36) 

The first term is an jS'-wave amplitude, the second is a longitudinal D-wave term, and the 
third is a P-wave term. We find for — > decay 

Ai = aifpmp{mD + mK*)Af^'' (ml) + a2fK*mK*{mD + mp)Ai^{mj^.), 

A-i = ^ aifpmpA2^* {ml) a2fK*mK-A2''{m\*), (37) 

rriD + rriK* mo + mp 

V = 1 ai/,m,V^^^*(mJ) a2fK'mK*V''f{ml.), 

rriD + rriK* mo + rUp 

where we have applied Eq.(29). Experimentally, a mixture of transverse and longitudinal 

polarization is found to be consistent with a pure 5*- wave term [22,18]: 

B{D+ K*\+- S - wave) = (4.11^;^) %, 

B{D+ K*\+; P - wave) < 5 x 10"^ (38) 
B{D+ K*\+] D - wave longitudinal) < 7 x 10^1 

It is convenient to write the Lorentz invariant amplitude in terms of three helicity amplitudes: 



^(m^ - m^, - ml)Ai + m\,plA2 
= A^ + rriDPcV, (39) 



H__ = Ai- tudPcV, 
where Pc is the cm. momentum. The decay rate is then given by 



r(D+ ^ i?*v+) = ^^^(i^oor + i^++r + \H-\'). (40) 



Assuming ri = r2 as before, Af^(0)/Af^*(0) = 0.89 [23], and a monopole behavior for form 
factors [m* = 2.53 GeV for (sc) current and 2.42 GeV for (uc) current], we obtain 

^{D ^ K*p) ^ -0.34:, riD ^ K*p) ^ -1.35. (41) 

The value of ^ is sensitive to the ratio Ai'^{0)/A^^* (0) which has not been measured. 
Nevertheless, it does not affect the general feature that \r{D ^ W)\ ~ \r{D ^ VP)\. 

To conclude this section, the parameter r2, which measures the nonfactorizable contribu- 
tion from color octet currents, is found to be of order ~ —2/3, — (0.9 ~ 1.1), — (1.2 ~ 1.3) 
respectively for D — > PP, VP, VV decays, in accordance with the expectation that soft- 
gluon effects become stronger when the final-state particles become less energetic. The 
parameters ai and especially a2 are thus not universal; they are channel or class dependent. 



3. Bottom Meson Decays In the exclusive two-body nonleptonic weak decays of bottom 
mesons, the parameter oi can be estimated directly from neutral B decays, e.g. B^ — > 
D+7r~, D^p' . This comes from the fact that, contrary to the charmed meson case, both 
W exchange and final-state interactions are presumably negligible as the decay particles are 
very energetic (for a recent study of final-state interactions in B decays, see Ref.[24]). As 
for the parameter 02, -B — > ^K^*^ are thus far the only color-suppressed modes which have 
been measured experimentally. Hence, a direct measurement of 02 is available only in these 
decay modes. A fit to the CLEO II data of B' 7pK-^*\ B^ jpK^^*^ yields [8] (see also 
Ref.[25]) 

|a2(^ ^ V'i?)! = 0.26±0.02. (42) 
Since the Wilson coefficients at the renormalization scale p ~ mb are given by 

ci(mfe) ~ 1.11 , C2(mb) ~ -0.26 , (43) 
it follows two possibilities: either 02 < and 

6(5 -^i^K)^0, r2{B^ ^pR) ^ -| (44) 
which is precisely what expected from the leading 1/A^c expansion, or 02 > and 

6(P ^ i^K) ^ 0.47 , r2{B i^K) ^ 0.27 . (45) 
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Unfortunately, what is the sign of 02 for B ipK^*^ decays is still very confusing. On the 
one hand, a recent calculation based on QCD sum rule indicates a destructive interference 
between the nonfactorizable soft-gluon term and the subleading 1/Nc factorizable term [13]. 
This in turn implies a negative 02. On the other hand, the ratio 02/01 is found to be positive 
in B- Dn, Dp, D*n, D*p decays [8]. Q 

In the following, we will argue that 02 appears to be channel dependent: Its sign is 
opposite in 5 — > -ipK^*^ and B D^*^n(p) decays. Following the argument presented in the 
last section, it is expected that 

\r2{B ^ ijK)\ > \r2{B ^ Dn)\. (46) 

A fit of the theoretical calculation [26] 

- BiBO D^n~) + ^-^^^J ^^^^ 

to the measured value i?i = 1.89 ± 0.26 ± 0.32 [8] gives Q 

— (E-^Dtt) ^ 0.30. (48) 

When combining with Oi ~ 1.05 determined from D~^7r~ , this leads to [] 

a2iB ^ Dn) ^0.32, ^(5 ^ Dtt) ^ 0.52 , r2(5 ^ Dtt) ^ 0.36 . (49) 

Evidently, the relation (46) can be satisfied only if 02 (-B — >• ipK) is negative: 

02(5 ^ ijK) = -(0.26 ± 0.02). (50) 

This argument that the nonfactorizable term and the subleading l/N^ factorizable term in 
B —>■ tpK decay are in opposite signs is also in accordance with a recent QCD sum rule 
calculation [13]. 

^The ratio 02/01 = 0.23 ± 0.11 obtained in Ref.[8] actually comes from (i) the determination of |ai| = 
1.15 ± 0.11 from B° D'-*^tt{p) decays, (ii) |a2| = 0.26 ± 0.02 from B ipK^*'^ decays, and (iii) the 
constructive interference in B- D'-*K{p) decays. A priori there is no reason to expect that \a2{B 
D(*)7r(p))| = \a2iB ipK^*'>)\ as the cm. momentum in B ijjK is 1683 MeV, while it is 2307 MeV in 
B Dn decay. 

direct fit of 02/01 to the other ratios R2, R3, R4 (for a definition, see Ref. [26]) also gives a value 

substantially larger than 0.23 (see Eqs.(27-30) of Ref.[8]). 

'^Unlike the charmed meson decay, r2 here is supposed not to vary significantly from B — > Dtt to B ^ D* p 

decays. 



11 



The most striking feature observed by CLEO [8] is that the parameter r2{B D^*^7c(p)) 
is positive, recalhng that the analogous quantity is always negative in charm decay. If our 
above argument is correct, we have to understand why the sign of r2 changes dramatically 
from B il)K to B ^ D^*^n(p) decays. This should be checked by lattice and QCD sum rule 
calculations. We see in charmed meson decay that the parameters ri and r2 are empirically 
close. Is this still true in bottom decay? The QCD sum rule calculation in Ref. [12] indicates 
ri{B^ D^Ti^) ^ with a sign opposite to r2{B^ Dn). However, a direct fit to the 
data of B'^ D^*^Tx{p) (except for B^ D^p^) gives a value of ai of order 1.01 ~ 1.05 (see 
Ref. [24] and Table XX of Ref. [8]) with the neglect of W exchange and final-state interactions. 
It thus appears that the sign of ri{B^ D~^tt~) could be the same as that of r2{B^ Dti). 
This issue should be clarified in the future. 

4. Conclusions We have considered the nonfactorizable contributions to exclusive two- 
body nonleptonic weak decays of charmed and bottom mesons arising from color octet 
currents. From the data we have determined the parameter r2, which measures the soft 

gluon effects, to be 0.67, -(0.9 ~ 0.11), -(1.2 ~ 1.3), ~ 0.36 respectively for D 

Kit, K*'k, K*p, and B Dtt decays. Therefore, the soft-gluon effects become more impor- 
tant when final-state particles become less energetic, as expected. In the l/Nc approach, the 
subleading contribution is characterized by the parameters = {l/Nc + ri/2) {i = 1, 2). It is 
evident that the coefficients ai (= ci + C,iC2) and especially 02 (= C2 + ^2Ci) are not universal; 
they are decay mode or class dependent. 

The very striking feature observed by CLEO that r2 is positive for B Dn decays, 
whereas it is always negative in charmed meson decays, needs to be checked by lattice and 
QCD sum rule calculations. We have argued that for B -ipK^*^ decays, the nonfactorizable 
term contributes destructively with the 1/Nc factorizable term so that the sign of 02 is 
negative. As a consequence, the leading 1/Nc expansion turns out working well for B — » 
ipK^*\ though it fails for B D^*^'k{p) decays. 

®This number is taken from Ref.[27], which is smaUer than the original estimate ri ^ —1.5 given in 
Ref. [12]. 
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